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CN ■ We compute a variety of two and three-point real-time correlation functions for 
>■ ' 

T-|- ' a strongly-coupled non-relativistic field theory. We focus on the theory conjectured 



to be dual to the Schrodinger-invariant gravitational spacetime introduced by Bal- 
asubramanian, McGreevy, and Son, but our methods apply to a large class of non- 
r^ I relativistic theories. At zero temperature, we obtain time-ordered, retarded, and 

^^ I Wightman non-relativistic correlators for scalar operators of arbitrary conformal 

dimension directly in field theory by applying a certain lightlike Fourier transform 
to relativistic conformal correlators, and we verify that non-relativistic AdS/CFT 
reproduces the results. We compute thermal two and three-point real-time corre- 
lators for scalar operators dual to scalar fields in the black hole background which 
is the finite temperature generalization of the Schrodinger spacetime. This is done 
by first identifying thermal real-time bulk-to-boundary propagators which, when 
combined with Veltman's circling rules, yield two and three-point correlators. The 
two-point correlators we obtain satisfy the Kallen-Lehmann relations. We also give 
retarded and time-ordered three-point correlators. 
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1 Introduction and Summary 

The AdS/CFT correspondence [Il[2l[3] has enjoyed much success in elucidating aspects of 
strongly-interacting relativistic field theories over the past twelve years. However, there 
exists a broad array of strongly-coupled systems which admit a field-theoretic description 
but do not exhibit Lorentz invariance. Cold atomic gases in the so-called unitarity limit 
provide a prime example of this. The unitarity limit refers to a regime in which the 
system resides precisely at the crossover point between a gas of Cooper pairs and a Bose- 
Einstein condensate. This regime arises when the interatomic potential is tuned such that 
the range of the potential vanishes while the scattering length diverges, thus removing all 
scales in the problem. The ability to realize this scaling limit experimentally by exploiting 
the phenomena of Feshbach resonances [1] has drawn wide attention to this system and 
emphasized the need for more progress on the theory front to balance the successes in 
the laboratory. Recent approaches to this problem that employ field theory techniques 
directly can be found in [5], [6l [71 18] . 

The unitary Fermi gas is a natural candidate for testing the waters of non-relativistic 
AdS/CFT as it possesses not just scale invariance but a larger group of symmetries known 
as the Schrodinger group, which can be thought of as the non-relativistic analog of the 
conformal group. The fact that the Schrodinger group is not only a close cousin of con- 
formal symmetry but that it can also be embedded in the conformal group was a key 
observation in formulating the first proposal for a non-relativistic AdS/CFT correspon- 
dence [HI [ID]- This proposal posited that the non-relativistic (Schrodinger) correspondence 
works in much the same way as its relativistic counterpart, but with a particular plane- 
wave-type geometry (sourced by additional non-trivial background fields) replacing AdS 
as the dual gravitational spacetime at zero temperature. 

The form of the Schrodinger spacetime suggests that we should think of the dual 
Schrodinger theory as being embedded within a parent theory living in one higher dimen- 
sion. If we denote the time coordinate in the parent theory by t, then the Schrodinger 
theory is obtained by switching to light cone coordinates x^ = y ±t, where y is one of 
the spatial directions, and then projecting the spectrum onto a sector of fixed x~ momen- 
tum. The x'^ coordinate plays the role of time in the Schrodinger theory. If x~ is also 
compactified, then the x~ momentum assumes a discrete set of values that are naturally 
interpreted as particle number for the non-relativistic Schrodinger theory, and the proce- 
dure is referred to as discrete light cone quantization (DLCQ). At zero temperature, there 
are subtleties involved with this compactification since x~ is lightlike, and one needs to 
consider large particle number in order to trust the gravity dual [TT| [T2| [13] . This idea 
of embedding the Schrodinger theory within a parent theory is closely related to the no- 
tion of Bargmann spacetimes [T^ [T5| [TB] , and it greatly facilitates the formulation of a 
non-relativistic version of AdS/CFT. 



The original conjecture of Refs. [HI [ID] has since been placed on firmer ground by 
subsequent work which uplifted the dual gravitational background to a solution of IIB 
supergravity and extended it to a finite temperature black hole solution [HI [12], [13] . Fur- 
thermore, it was shown that at finite temperature, x~ is no longer a lightlike direction, so 
that subtleties pertaining to its compactification no longer arise. The work of [HI [121 [13] 
also clarified that the parent theory related to the Schrodinger metric constructed in [91 [TU] 
is an example of a well known class of non-commutative theories known as dipole theories 
[T71 [T8| [T9] . The dipole theories are obtained by starting from an ordinary commutative 
field theory and replacing regular multiplication of fields with a non-commutative star 
product. The particular parent theory described by the Schrodinger metric is the dipole 
theory that results when A/'=4 super Yang-Mills is deformed in this way, and the corre- 
sponding Schrodinger theory is just the DLCQ of this. In what follows, we will always 
refer to these theories simply as the parent and Schrodinger theorie^j, and we will also use 
these terms to distinguish between their gravity duals, which are also related by DLCQ. 
These discoveries constitute an important milestone toward extending AdS/CFT to the 
realm of non-relativistic theories and applying it to improve our understanding of systems 
like the unitary Fermi gas. 

Some of the most important information that AdS/CFT provides about a field theory 
are its correlation functions. As for theories exhibiting conformal invariance, the two and 
three-point correlation functions for Schrodinger-invariant theories also provide crucial 
checks of the conjectured correspondence since these correlators are largely determined by 
the symmetry group alone [2CT| [51] . In particular, the two-point functions are completely 
fixed by the symmetry up to an overall constant. Unlike the conformally-invariant case, 
however, it is not true that the functional form of the three-point functions are fully 
determined. It was verified in [22 l j23 | [23] that AdS/CFT reproduces both two and three- 
point scalar correlation functionso 

In this paper, we will focus on computing real-time correlation functions for scalar fields 
in Schrodinger field theories. Real-time correlation functions are of particular interest 
since these in turn yield quantities like conductivity and viscosity. However, real-time 
correlation functions in AdS/CFT pose more technical challenges relative to Euclidean 
ones since it is less clear how to formulate basic AdS/CFT recipes in Minkowski signature, 
where the bulk spacetime tends to be more complicated. Early attempts to deal with 
this problem yielded a case- by-case treatment [2S1 [2S], and only in the last two years 



^Although Schrodinger invariance is broken at finite temperature, we will continue to use the term 
"Schrodinger theory" even in this case. 

^Intriguingly, the authors of [35] also computed the three-point function using a certain field theory 

model for a unitary Fermi gas and showed that, in addition, AdS/CFT correctly computes the piece of 
the three-point function which is not determined by Schrodinger symmetry. 



have more systematic approaches appeared in the hterature [2Zl ESI [22] • In the context of 
Schrodinger-invariant theories, the methods of [271 ESj were employed by Hoang and Leigh 
[21] to compute time-ordered and Wightman two-point functions at zero temperature. We 
will show that the approach given in [29] reproduces their findings in a simple way, and 
we further use it to compute a wide variety of real-time two and three-point correlators 
at both zero and finite temperature. Our results for the zero-temperature two-point 
functions are also consistent with the bulk-to-bulk correlators computed in [3D] . 

The complete functional form of the zero-temperature two and three-point correla- 
tion functions can in fact be computed in real time without having to invoke AdS/CFT 
or any other method for solving the Schrodinger field theory due to the intimate rela- 
tion between the Schrodinger and conformal groups. More specifically, we may obtain 
correlation functions of the Schrodinger theory directly from the correlation functions of 
a CFT by switching to light cone coordinates and Fourier-transforming with respect to 
x~, as was originally noticed some time ago in Ref. [21] and more recently exploited in 
the context of AdS/CFT in Refs. [221 123]. This procedure is essentially equivalent to 
performing DLCQ, with the role of the Fourier transform being to project onto a sector 
of fixed x~ momentum. However, as stressed in [21], it is important to note that the 
Fourier-transform trick must be performed keeping x~ noncompact so as to avoid having 
to construct CFT correlators on a lightlike circle. We will still refer to this procedure as 
DLCQ since we are free to consider fixed particle number in the resulting expressions, 
keeping in mind the zero-temperature subtleties mentioned earlier. 

We will apply the Fourier-transform technique to compute various zero-temperature 
correlators in real time. Although the authors of [2^ [22] computed real-time correlators, 
they did not keep track of the type of correlator (e.g. time-ordered, retarded, etc.), as their 
primary focus was on testing whether non-relativistic AdS/CFT reproduces Schrodinger- 
invariant expressions. Since the ie prescriptions which distinguish between the different 
types of real-time correlators are well known in the case of a relativistic CFT (see [2H] for a 
review), it is a straightforward task to perform DLCQ on these correlators to produce real- 
time Schrodinger correlators, and we will show that the type of correlator is preserved 
under DLCQ. In addition, we verify that the standard Kallen-Lehmann relations are 
satisfied by the various real-time Schrodinger correlators. 

We stress that the CFT correlators that we Fourier transform are not the correlators 
of the parent theory. In particular, note that the parent dipole theory is not Lorentz- 
invariant [T71 [IS, [H] • Therefore, the resulting Schrodinger-invariant correlators will not 
contain the same overall constants as those of the DLCQ of the dipole theory. Our aim in 
applying the Fourier transform to CFT correlators is to ascertain the form of the various 
real-time Schrodinger correlators at zero temperature; this provides us with an important 
check of the exact Schrodinger theory correlators we will compute from non-relativistic 
AdS/CFT. 



We compute the zero-temperature correlators in momentum space as well as position 
space. DLCQ applied to momentum space correlators simply amounts to performing a ro- 
tation in the plane spanned by u and py, the energy and momentum associated with t and 
y. In the case of three-point functions, we express the time-ordered and retarded correla- 
tors in a way that is very reminiscent of AdS/CFT, namely in terms of an integral over 
what are readily interpreted as three bulk-to-boundary propagators in the Schrodinger 
spacetime. This allows us to define in a natural way non-relativistic Feynman and re- 
tarded bulk-to-boundary propagators which in turn are easily verified to reproduce the 
corresponding two-point correlators that we compute directly in field theory. Wightman 
propagators are then defined in such a way that they reproduce Wightman two-point 
functions. The exercise of constructing the different real-time bulk-to-boundary propaga- 
tors and establishing their interrelations not only provides additional consistency checks 
of non-relativistic AdS/CFT, but also helps set the stage for the computations of thermal 
correlation functions, which rely heavily on correctly identifying these propagators. 

One of the main objectives of this work is to compute real-time two-point and three- 
point non-relativistic correlation functions for scalar operators dual to minimally coupled 
scalars at finite temperature. The Schrodinger theory at finite temperature also possesses 
a nonzero chemical potential which breaks Schrodinger invariance, meaning that we can 
no longer receive guidance by applying the Fourier-transform trick to CFT correlators as 
in the zero-temperature case. With the confidence gained from cross-checking AdS/CFT 
with field theory at zero temperature, we therefore focus solely on AdS/CFT in the finite 
temperature case. 

The first step is to compute real-time bulk-to-boundary propagators. This is done 
by starting with the Euclidean version of the metric which describes a black hole in 
the Schrodinger spacetime [HI [121 [13] ^^^ solving the scalar wave equation to obtain 
the Euclidean bulk-to-boundary propagator in momentum space. This is the Euclidean 
propagator in the gravity dual of the parent dipole theory. We analytically continue 
this to a retarded propagator and then perform DLCQ to get the retarded Schrodinger 
bulk-to-boundary propagator. In [29], it was shown that in the relativistic case, the 
other real-time thermal propagators can be derived from the retarded one using certain 
identities. These identities naturally extend to the parent dipole theory, where we can 
then apply DLCQ to show that similar relations between Schrodinger propagators also 
hold. 

Once we have the real-time Schrodinger bulk-to-boundary propagators, we can assem- 
ble them into real-time two and three-point correlators. This problem was solved in the 
case of thermal A/'=4 SYM in [29] by extending Veltman's circling rules [311 [321 [33] to the 
gravitational theory on the black hole spacetime, and similar arguments should apply in 
the present context as well. We compute two-point correlators in this way and confirm 
that the results satisfy standard Kallen-Lehmann relations. The retarded three-point was 



worked out explicitly in [22], and we borrow the result to immediately write down the 
three-point function for the parent dipole theory. A simple application of DLCQ to the 
result in turn gives the retarded Schrodinger three-point correlator. In addition, we apply 
the circling rules to construct the time-ordered Schrodinger three-point function. 

Before proceeding with the calculations, we pause for a moment to briefly discuss the 
dimensionality of the theories we consider. We keep the number of spatial dimensions 
(denoted d) arbitrary as much as possible. There is certainly no obstacle to doing so in 
our discussion of DLCQ applied directly to CFT correlators at zero temperature. Further- 
more, we may still consider the d+3 dimensional Schrodinger spacetime as a "bottom-up" 
holographic model even though this spacetime can only be uplifted to a solution of JIB 
supergravity in the case d = 2. Leaving the dimension arbitrary here also facilitates com- 
parison with our field theory results. We could adopt a similar philosophy toward the 
Schrodinger black hole spacetime, however it is considerably more difficult to solve the 
wave equation for values of d not equal to 2, so we restrict ourselves to (i = 2 in this case. 

We will consider various types of correlators and bulk-to-boundary propagators through- 
out the paper, and we distinguish between these with different fonts of the letter g, as 
we will now clarify. Relativistic CFT correlators are denoted by g, while non-relativistic 
Schrodinger correlators are denoted by G. The symbol g specifies bulk-to-boundary prop- 
agators in the gravity dual of the parent theory. We also define bulk-to-boundary prop- 
agators for the gravity dual of the Schrodinger theory, and these we denote by ^. In 
addition, we use tildes to denote momentum space functions, the superscript (3) to de- 
note causal three-point correlators, and additional decorations to distinguish between 
the different types of real-time functions. In particular, the subscripts F, R, and A de- 
note time-ordered, retarded, and advanced functions respectively. Reverse-time-ordered 
functions are distinguished from time-ordered functions by an additional bar (e.g. Gp) 
which should not be confused with complex conjugation, for which we use the symbol *. 
Wightman two-point functions have superscripts ±, while Wightman three-point functions 
are specified by a subscript which shows the order of operator insertions. For example, 
6*123 is a Schrodinger Wightman three-point function in momentum space proportional 
to (O1O2O3), where the Oi are arbitrary scalar operators in the Schrodinger theory. 

The paper is organized as follows. In section [21 we compute all types of scalar real-time 
Schrodinger two-point correlators at zero temperature in both position and momentum 
space by applying the Fourier-transform trick to relativistic CFT correlators. We also 
verify that the results satisfy the Kallen-Lehmann relations. In section [31 we perform 
similar computations for real-time zero-temperature three-point correlators. We confirm 
that standard relations among the three-point functions hold, and we also express the 
results in a way that is very suggestive of an AdS/CFT calculation, allowing us to identify 
Schrodinger "bulk-to-boundary propagators" . In section [H we construct real-time bulk- 
to-boundary propagators directly in the zero-temperature gravity dual and show that 



the results are consistent with the answers obtained using the Fourier-transform trick 
in sections [2] and [31 Finally in section we apply the intuition developed in the zero- 
temperature case to construct thermal real-time bulk-to-boundary propagators for scalar 
fields. These are then employed to compute two and three-point scalar correlators for the 
Schrodinger theory at finite temperature. The first four appendices each contain details 
about the application of the Fourier-transform trick to particular zero-temperature three- 
point correlators considered in section [31 In appendix [Ej we construct an explicit example 
of a massive scalar in the effective 5d gravity theory which arises when a fluctuation in a 
particular off-diagonal metric component has non-trivial charges on the 5d internal space. 



2 Zero-temperature two-point functions from CFT 
correlators 

In this section, we will obtain Schrodinger-invariant non-relativistic zero-temperature two- 
point functions in position space by starting with real-time relativistic conformal two-point 
functions in light cone coordinates and Fourier-transforming with respect to one of these 
coordinates. We also obtain momentum space two-point functions by translating the 
Fourier-transform trick to momentum space, where it becomes a simple redefinition of 
momenta. We will pay particular attention to how the different types of relativistic 
real-time functions (e.g. time-ordered, retarded, etc.) map to the different types of non- 
relativistic functions. We will see that the type of real-time correlator is preserved under 
the special Fourier transform. We also verify that standard Kallen-Lehmann relations 
between the various real-time functions are satisfied by the non-relativistic correlators we 
obtain. 

2.1 Time-ordered and reverse-time-ordered two-point functions 
in position space 

Consider first the time-ordered relativistic conformal function: 

gF{x^,x) = -—T—— 5 — X- (2-1) 

Here, x^ = y ±t, where t is the relativistic time coordinate, and y is one of the spatial 
coordinates. We use the mostly plus signature for Minkowski space. The coordinate 
vector X represents the remaining d spatial coordinates which coincide with the spatial 
coordinates of the non-relativistic theory: x = {xi,X2, ■■■,Xd). Following [211 [23], we can 



obtain a non-relativist ic two-point function Gp by performing a Fourier-transform with 
respect to the x~ direction: 

This integral can be done by introducing a Schwinger parameter: 

Gp = -^^^^ r dss""-' I rfa;-e-[-^(-+--+-')]^e-^^^" . (2.3) 

The integral on x~ yields a ^-function: 



Gf = V-FTaT / dss^-'^-^'-" ^'^iM - x^s). (2.4) 

27ri (A) Jo 



The 5-function will evaluate to zero unless Mx"*" > 0. Therefore, we find 



-t 



vA+l iM{x^+ie)/x-* 



Gf{x\x) = ^-^|Mr-^g(Mx+) ^^^^^ ■ (2.5) 

This has the form of a non-relativistic time-ordered two-point function where the non- 
relativistic time coordinate is x^ /2. Hoang and Leigh [21] found the same result employing 
a very different approach [271 [28]. We have retained the ie in the final expression because 
it removes a potential singularity at a;"*" = 0. It is clear from (12. 2p that this singularity 
should not be present. 

If we instead start with the reverse-time-ordered relativistic two-point function: 

—i 1 
gF{x^,x) = -g*p = ^. — ^ r^^, (2.6) 

the same calculation now gives 






^^^^ '^^ = ^^^f(A)'^' ^^"^^^ |^+|A ■ (2-7) 

This has the form of a non-relativistic reverse-time-ordered function. 



2.2 Time- ordered and reverse-time-ordered two-point functions 

in momentum space 

To make it clear how the Fourier-transform trick translates to momentum space, we first 
consider a simple example. The Fourier transform of the conformal time-ordered two-point 

8 



function with A = 1 and d = 2 is 

~gp{u^,Py,p) = f^ j dtdyd'xe^-'-^^yy^^^^-' / = / „ . , (2.8) 

471"^ y — t"' + y^ + x^ + le u^ — Py — p^ + le 

where p = {pi,P2) and p = \p\. Introducing the hght cone coordinates x^, we can write 
this as 



9f 



dx~^ dx~ d'^xe^^^'^~^y^^'^^^^ '^^^'^^''^y''^'^^^ -jp-x_ 



Stt^ J x^x + x'^ + ie 

We now define the non-relativistic energy Q and mass M in terms of the relativistic energy 
u and momentum component py-. 

^^liUJ-Py), M^^{u+Py). (2.10) 

Applying the inverse of this map, 

u; = M + n, py = M-n, (2.11) 



we obtain 






1 1 



2M2(]-|i + ,eM 



(2.12) 



This is the non-relativistic time-ordered propagator in momentum space with non-relativistic 
energy 

-E'non-rcl. = 2^2. (2.13) 

It is easy to check that (I2.12p is the Fourier transform of (12.51) with A = 1 and d = 2jj 

Gpin.p) = i I dx+d'^xe'^^^-'^-^Gpix^^x). (2.14) 

Obviously the same approach can be applied to any other momentum-space function 
since, from the point of view of the Fourier transform in (12. 8 p and (12.90 . we are simply 



■^Our conventions for both relativistic and non-relativistic Fourier transforms are established in this 
subsection. In particular, note that we define our non-relativistic Fourier transform with an additional 
overall factor of i because the non-relativistic time coordinate is .t+/2. 

9 



redefining momentum-space variables. Starting from tlie relativistic time-ordered two- 
point function witfi arbitrary A and d, 

9F[UJ,Py,P) 



47^2 



dtdyd xe 



2™„«(^t-ipH?/-ip-^. 



(-t2 + 2/2 + 2;2 + ze)^ 



(2i)'^-2^7ri-ir(f + 1 - A) 
f(A^ 



/ 2 2 2 , ■ \ 

[u -Py-p +te) 



A-l-1 



(2.15) 



the map (12. lip gives the non-relativistic time-ordered functioiu: 

~. ^ (2z)'^-2A^i-ir(| + 1 - A) , , ,A-^-l 



r(A) 



^ ^^ -\M\^^2-^ 

2^-i+ir(A) 



2fi 



2M 



+ ieM 



A-f-l 



(2.16) 



The reverse-time-ordered two-point function is just minus the complex conjugate of Gp'- 

^F{n,p) = -G*p 



g-m[(f-A)e(M)-e(-M)]^|-ipCd + 1 _ A) 

-, -\M\ 

2^-|+ir(A) 



A-#-l 



2fi 



J9^ 

2M 



ieAf 



(2.17) 



2.3 Wightman two-point functions in position space 

Next consider the following relativistic Wightman function: 



g^{x^,x) 



1 



Air"^ (x+x~ + x'^ + ie(x+ — x~))^ 
We can facilitate taking the Fourier transform by rewriting this as 

-. -^ 1 

9 



(2.18) 



(2.19) 



Att"^ {{x+ — ie)x~ + x"^ + ie{x+ — ie))^ 

Before taking the Fourier transform, it is useful to first factor out a factor of x~^ — ie from 
the denominator: 



GHx\x) = ^{ 



x~^ — ie) I dx 



-iMx 



ix- + -^ + le)^ 



(2.20) 



'*The second line of ([2l6| can be obtained with the help of the identity (a6)" = Qa^ag27rja,,(a,6)^ 
7?(a,&) = 6l(-Im(a))0(-Im(6))6'(Ini(afe)) - 6'(Im(a))6l(Ini(fe))e'(-Ini(afo)). 
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One can then use the positivity of e to introduce a Schwinger parameter as we did in the 
case of time-ordered functions in the previous subsection. The x~ integral then gives a 
(5-function: 

G+ = ^"'^^^l-rfAV^'^ r ^^^^"'e-[^-*^'/(^^-^^)]^5(M - s). (2.21) 

Since the argument of the 5- function does not contain x^, a step function involving x^ 
does not arise in this case. The result is then 

This agrees with the result found in [23] using the methods of [271 ESj • The role of the ie 
in the exponent is to ensure that G^ vanishes at x"*" = 0, while the ie in the denominator 
controls the location of the branch cut for non-integer A^ 

(x+ - ze)-^ = e{x+)\x+\-^ + e(-x+)e*"^|x+|-^. (2.23) 

If we instead consider starting with the other relativistic Wightman function, 

g-{x^,x) = ZL ^ (2.24) 

47r^ (x+x -|- x^ — «e(x+ — X j)^ 

the same steps lead to an s-integration as above, but with biM + s) instead of (5(M — s) 
appearing in the integrand, yielding 

■A-l iMx2/(x++ie) 

Using (12.51) . (I2.22p . and (12.251) it is easy to check that the relations 

Gf[x^,x) = e{x+)G+{x+,x) + e{-x+)G-{x+,x), 

Gf{x+,x) = e{-x+)G+{x+,x) + e{x+)G-{x+,x), (2.26) 

are satisfied|j 

^Throughout this paper, we adopt the convention that (—1) = e*'^ if no ie prescription is given 
exphcitly. 

^When x+ 7^ 0, the ie in the exponents can be discarded relative to a;+. When a;+ = 0, Gp, G+, and 
G~ each vanish identicahy, and p.26p is satisfied triviahy. 
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2.4 Wightman two-point functions in momentum space 

Now consider the relativistic Wightman functions with A = 1 and d = 2 in momentum 
space: 

~g^{co,Py,p) = -2m9i±tu)Situ''-pl-p''). (2.27) 

Under (12.111) . these transform to 

G^{n,p) = -2nie{±{n + M))6{AMn - p"^). (2.28) 

For M > (M < 0), the (5-function has support at positive (negative) values of Q, so that 

G+(fi,p) = -2me{M)5{AMn-p'^), (2.29) 

which is the Fourier transform of 02.221) for A = 1, d = 2, and 

G^{n,p} = -2'Kie{-M)5{AMVL - p^). (2.30) 

To obtain the non-relativistic Wightman functions for general A and d, we start with 
the Fourier transform of (I2.18P : 

The map (12. lip then gives 



2rf+l-2A^- 



G^{^,P) = - V(A)r(A - i^ (^(M)e{AMn-p'){AMn-pY~-^~'- (2.32) 

We have used that 

0{co - ^Jpl+p') = eicu + py)e{cu' -pI- p') (2.33) 

to rewrite the argument of the step function before applying the map. The second rela- 
tivistic Wightman function has the form 

2rf+l-2A^| 



r{co,Py,p} = -' r(A)r(-A-'-) ^^~"' ~ \/p'+p')(^'-pS-p')^"^"'- (2-34) 



r(A)r(-A-f 

Now observing that 



0{-io - ^pl+p') = e{uj' -pI- p^i-uj - py), (2.35) 
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we see that the map (12.111) leads to 



d 



G-{n,P) = - V(A)r(A - ^) ^(-^)^(4^^ - p')(^Mn - pY''^^'- (2.36) 

So far, we have seen that the relativistic (reverse-) time-ordered two-point function trans- 
forms into a non-relativistic (reverse-) time-ordered two-point function, and the relativistic 
Wightman functions map to non-relativistic Wightman functions. It remains for us to 
check what happens to relativistic retarded and advanced functions under this mapping. 

2.5 Retarded and advanced two-point functions in position space 

The relativistic conformal retarded two-point correlator can be written in terms of Wight- 
man functions: 

gnit, y, x) = e{t) [g^{t, y, x) - g-{t, y, x)] . (2.37) 

In terms of light cone coordinates, we therefore have 
9r{x^,x) = —e{x^-x~ 



((x+ — ie)x~ + x'^ + ie{x+ — ie)))^ 

{{x+ + ie)x- -h x2 - ie{x+ +ie))^\' 

Factoring out the (x"*"=F2e) in the denominators, assuming x"^ 7^ 0, and Fourier-transforming 
leads to 

—i f^^ ^-iMx- 

Gr(x^,x) = -—r{x^-ie)-^ dx' 



A-k"^ 7_oo {x~ + x"^ / x^ + ie)^ 



i , J. . sA 



-iMx- 



+ --(x+ + ie)^/ dx~- ^-- —J-. (2.39) 

Notice that we have chosen to write x^ jx^ instead of x^/(x"^ ± ie) in the denominators. 
Since we are assuming x'^ 7^ 0, we may rewrite x^/(x"*" ± ie) = x"^ /x^ =F iex"^ / {x^Y , 
and the term iex"^ /{x'^Y can be eliminated by rescaling the e in the denominators of the 
integrands. The ie in the factors (x^ ± ie)~^ appearing outside the integrals cannot be 
neglected even though x"*" 7^ because the ie in these factors control branch cuts, as we 
discussed in the previous subsection on Wightman functions. 
Changing the integration variable to u = x~ -|- x^/x"*", we have 

_■ i>x'^/x++x+ -iMu 

Gj, = _4(x+-ze)"V^^'/^W du- 



4vr2' i-oo (n + ze)^ 

i-x'^/x++x+ „-iMu 



+ -^(x+ + ^e)-^/ du- —^. (2.40) 

47r2^ i-oo iu-ie)^ 
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The integrals can be evaluated using 



and we find 



!Toy_''"FTi;F = ±^i^A'"-'''(")''(±A^). (2-«; 



Gb{x*,x) = i 'l'^l\ |M|^-'e"'^-'""-"'0(ir+)|ir+|- V'"'-'°'% 1+ ^ 0. (2.42) 

27rr(A) 

This has the form of a non-relativistic retarded two-point function. This expression is 
valid for x'^ 7^ 0; when x'^ = 0, it is easy to check that G^i = 0. (This can be done by 
setting x"*" = in fl2.38p and Fourier-transforming.) An analogous calculation reveals that 
if we start with the relativistic advanced function, we obtain the non-relativistic advanced 
function: 

Ga{x+,x) = ^ X^ J M|^-^e^"(^-^)'(^^)g(-x+)|x+|-V^"'/"^ x+ ^ 0. (2.43) 
27rr(A) 

A shortcut to obtaining this result is to write down the analog of 02.391) for Ga and to 
notice that Ga can be obtained from G^ by sending x~^ — )■ —x~^ and M — )■ — M in the 
final answer. It is easy to check using (12.221) and (I2.25P that standard relations between 
retarded/advanced and Wightman functions: 

Git = e{x+)[G+ -G-], GA = 0i-x+)[G- -G+], (2.44) 

hold for the non-relativistic two-point functions we have obtained. 

We conclude that the type of real-time two-point functions is preserved under the 
mapping which takes us from relativistic to non-relativistic functions. 



2.6 Retarded and advanced two-point functions in momentum 

space 

The Fourier transform of the relativistic retarded two-point function (I2.38P with A = 1 
and d = 2 is given by 

gR{u,Py,p} = 9{u)gF + e{-u)~g*p = — . . (2.45) 

u^ — p^ — P'^ + leuj 

Applying ( 12. lip , we obtain 



AMVt -p^ + ie AMVt - p^ - ie 
^ (2.46) 



AMVL - p"^ + ie{VL + M) 
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This function has a pole at 



"' = £j-ni+liM5l- (247) 



Since this pole lies in the lower half fi-plane regardless of the value of M ot p^, we see 
that this function has the analytic properties expected of a retarded Green's function. 
We may therefore redefine e to absorb the positive coefficient which multiplies it in Qpi 

GR{n,p) = —, . (2.48) 

The advanced function is the conjugate of this: 

Ga = G*j,. (2.49) 

Comparing these results with (I2.12p . we see that Gp, Gr and Ga obey the following 
relation: 

Gpisi.p) = e{M)GRin,p} + e{-M)GAin,p}. (2.50) 

On the other hand, we can obtain what appears to be a slightly different relation directly 
from (EJSD: 

Gpi^,?} = o{n + M)GR{n,p} + e{-n - m)Ga{^,p). (2.51) 

At first glance, it may seem surprising that these two relations hold simultaneously. 
However, their compatibility has a simple origin in the form of the imaginary part of 
the retarded two-point correlator. Using (12.461) . we can write this as \yivGr ~ sgn(r2 + 
M)(5(4Mr2 — p^). Note that only the imaginary part is affected by the presence of the 
step functions in (I2.50p and (12.511) . Since the (5- function requires that Vt and M have the 
same sign, this is equivalent to IovGr ~ sgn(M)(5(4Mr2 —p^), which agrees with (12.481) . 
Also notice that we can write the relation in a third way: 

Gf{q,p) = 9{n)GRin,p) + e{-n)GA{n,p). (2.52) 

This form of the relation makes it apparent that there is no chemical potential at zero 
temperature. We will now see that the relations (12.501) . (I2.5ip . and (12.521) hold for general 
values of A and d. 

The Fourier transform of (12.381) for general A and d is 

gR{u, py, p) = 9{uj)gF + 6{-u)g*p, (2.53) 
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where gp was given in fl2.15p . This yields 



GRi^.p) 



9{n + M)Gf + 9i-n - M)G*p 
(2zsgn(fi + M))'^-2^7ri-ir(f + 1 - A) 



r(A) 



{AMVt - p'^ + ieiyt + M))' 



£'"1(2 



i7r[(f-A)e(A/)-e{-A./)] 



Trf-irr^ 



1-A) 



2A- 



-ir(A) 



|M|^- 



2fi- 



P 

2M 



+ ie 



A-l-i 



(2.54) 



Similar reasoning leads to the following advanced two-point function: 



Ga{^,P) 



9{n + M)G*p + 9{-n - M)Gf 
i-2isgn{n + M))'^-2^7ri-ir(f + 1 - A) 



r(A) 



{AMQ - p"^ - ie{n + M)) 



A-f-l 



-J7r[(|-A)0{J\/)-6»(-J\/)] 



~l-pfd 



712 T^ + 1 



A) 



+1 



r(A) 



\M\ 



A-f-l 



2f] 



2M 



2e 



A-l-i 



(2.55) 



It is not difficult to check that the relations (I2.50p . (I2.5ip . and (12.520 hold for general A 
and d. Furthermore, the relation (I2.50p provides a useful check of the consistency of our 
results since it must also hold in position space, and one can readily verify that ( 12.50 . 
(I2.42p . and (12.431) satisfy (I2.50p . We will see later on in section [5] that the second relation, 
(I2.5ip . generalizes to a relation that is also valid at finite temperature. 



3 Zero-temperature three-point functions from CFT 
correlators 

In this section, we apply the Fourier-transform trick to relativistic three-point functions 
to obtain Schrodinger-invariant non-relativistic three-point functions in position space. 
We do this for time-ordered and Wightman functions and show that the usual identity 
relating these types of correlators is satisfied. In position space, retarded and advanced 
functions can be expressed in terms of the Wightman functions [36] . 

We compute non-relativistic time-ordered and retarded/advanced three-point func- 
tions in momentum space as well. This can be done either by computing the Fourier 
transforms of the non-relativistic position space three-point functions, or by starting with 



16 



relativistic conformal three-point functions in momentum space and performing an ap- 
propriate redefinition of momenta as we did for two-point functions in section HJ We 
will apply the former approach to obtain time-ordered functions and the latter approach 
for retarded/advanced functions. Wightman three-point functions are quite unwieldy in 
momentum space, so we do not include them here. 

The final answers we obtain for the momentum space three-point functions are of 
the form one would expect from an AdS/CFT calculation. In particular, we express the 
results as integrals of products of three functions which are naturally interpreted as bulk- 
to-boundary propagators for scalar fields in non-relativistic AdS/CFT [9l[T0]. In section 
m we reproduce these bulk-to-boundary propagators by applying DLCQ to propagators in 
the gravity dual of the dipole theory and show that they give rise to two-point functions 
which are consistent with the non-relativistic correlators we computed in section [21 

3.1 Time-ordered three-point functions in position space 

The relativistic time-ordered three-point function in position space has the form 

9^? = -{0\TO,{x,)O,{x,)Os{xsm = - . + - .\ .^.s. 

["''23"''23 "•" -^23 "•" ''^J 

(3,1) 



Fl3-^13 + ^13 + ^^J ^ Fl2^12 + ^12 + ^N 

We have introduced the notation Xij = Xi — Xj. The constants 6i are related to the 
conformal dimensions Aj of the scalar fields: 

<5i = ^(A2 + A3 - Ai), 62 = ^(Ai + A3 - A2), 63 = i(Ai + A2 - A3). (3.2) 

(3) 
To obtain the non-relativistic correlator Gp , we must Fourier-transform with respect to 

the x~ . The integrals can be computed with the help of Schwinger parameters, and the 

details can be found in appendix |Xl The result is 

Gf = -M^til!^^:^^ 

(3.3) 
where 

9 . o . 9 . 

X19 ~r t-c ^ort \ t-c -^i Q ~r t-c y < 

^^12 = ''+ + ^ -^^^ (3-4) 

•^12 "''23 "''13 
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and 



/ = e{xi,)e{xtM4,)h++ + ^(x+ )^(x+3)^(x3+ )/++- + e{xtM4iW{xt2)h- 

+ ^(X+ )^(X+ )^(X+3)J_++ + ^(X+ )^(X3+ )^(x+3 )/__+ + ^(X+)0(X+)^(X+ )/___. 

(3.5) 
The functions lijk appearing in this expression are given by 

+ e{-M2)e{Mi - |M2|)(Mi - iMsD'^^+'^^^^lMal'^^-^e^i^l^-^^l 

X B{5^,52)^i{5^A - 5,,5^ + 82 }^^}' ^\ iv^2{M,-\M2\))\ (3.6) 

\M2\ ) 



poo 
l"''12l 



/■oo 

/ dzz^'-^\M2 + ;z|'^^"^|Mi - zl^^-^e'^'^'OiMi - z)e{-M2 - z) 
Jo 

xt2\~^'0{Mi)e{-M2)\e{Mi - \M2\)M^^'-'\M2\^'+^'-'B{63,6i) 



X $1(53, 1 - 62, <5i + 53, ^, tVu\M2\) + e{\M2\ - M,)M','+'-^~'\M2\''~' 

X B{5s,S2)MS3,l-Si,52 + 5s,^,tvi2MA, (3.7) 



/■CO 

x^al"^' / ^^^^'"^2 - z\^'-^\Mi + zl^'-^e-'^'^'eiMi + z)e{M2 - z) 
Jo 

xt2\~^'9{-M2)\9{Mi)9{\M2\ -Mi)e^"^2^^iMf3-^(|M2| - MO'^^+'^^-i 



X B{62, <5i)$i(<52, 1 - ^3, 6^ + ^2, ^^^^^, ^v^2{\M2\ - M,)) 

+ e{-M^)\M^\'^-^\M2\''^''''B{53,5^)^^{53A - 82,5^ + 5-i,-^^,iv^2\M2\)\ 

(3.8) 

The function B{x,y) is the Euler beta function, while $1 is a hypergeometric function 
of two variables. These functions are defined in ( IA.13I) and ( lA.M^ respectively. The 



remaining three lijk appearing in (13. 5 p can be obtained from the ones quoted above using 



the relations 

J_++(5i,52,(53,Mi,M2,Wl2) = /+++(52, 51,^3,^2, Ml, -t;i2), 

/__+(5i, 52,^3, Ml, M2,i;i2) = /++_(<52,5i,<53,M2,Mi,-t;i2), 
/___(5i,52,53,Mi,M2,t;i2) = /+— (52,'5i,^3,M2,Mi,-t;i2). (3.9) 

We will see in section 13.31 that the lijk are proportional to the Wightman functions Gijk 
in precisely the way one would expect from the standard identity 

Gf = e{xt2)e{xtM43)Gi23 + ^(42)^(43)^(42)^132 + ^(42)^(4)^(42)^312 

+ ^(4)^(4)^(43)^213 + ^(4)^(4)^(43)^231 + ^(4)^(4)^(42)^321. 

(3.10) 

3.2 Time-ordered three-point functions in momentum space 

We next consider time-ordered three-point functions in momentum space. We could 
compute these by starting with conformally-invariant time-ordered three-point functions 
in momentum space and performing a redefinition of momenta, but we choose instead 
to Fourier-transform the CFT correlator in position space, eqn. 03.11) . As in previous 
sections, the Fourier transform is computed by introducing Schwinger parameters; the 
details can be found in appendix [Bl The result is 

X f dzzi JJ(-4M,fi, + p] - leY^'^K,^ (^z^ -AM^Vt^ + p] - le^ , (3.11) 

where we have defined 

z.,^5^<5,-^-l = A,-^-l. (3.12) 

Our result for Gp has precisely the structure one would expect to see arise out of an 
AdS/CFT calculation. This observation suggests that we can identify a non-relativistic 
time-ordered "bulk-to-boundary propagator" : 

Qpiyt^V, M, A) ~ zi+i(-4Mfi +p2 _ zey/^K^ fz^J-AMn+p^-ie) . (3.13) 

z can be interpreted as a radial coordinate in (a modified version of) AdS. We have put 
in by hand the factor of 22+^ appearing in Qp; this will be justified later by invoking the 
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proposed AdS/CFT dual O [ID], which we will review and make use of in section HI Note 
that after separating out a factor of z2+i for each propagator, there remains the factor 
z~'^~^, which is precisely the measure factor of the dual gravity spacetime. 

We can work out the position-space analog of ( 13.111) by computing its Fourier trans- 
form, with the result 

X Qf{Xi - X, Z\ Ml, Ai)^ir(x2 - X, Z\ Mg, A2)^f(x3 - X, Z\ M3, A3), 

(3.14) 

where 

gFix,z;M,A) = ^ X ^ , |M|^"^^^g(Mx+)|x+|-^e'*^("'+"')/"^ (3.15) 

27rr(A) 

can be interpreted as a Schrodinger Feynman bulk-to-boundary propagator in position 
space in the sense that it is the Fourier transform of fl3.13p . A similar expression was 
obtained by Volovich and Wen [23] . 



3.3 Wightman three-point functions in position space 

A conformally-invariant Wightman function in position space can be expressed in the 
form 

gi23{xi,X2,X3) = -(0|Oi(Xi)O2(x2)O3(x3)|0) = 



^23^23 + ^23 + ^^(,^23 ^23)J ^ 



Fl3^13 + ^13 + *^(^13 ^13)] " Fl2^12 + ^12 + ^^(^12 ^12)] ^ 

(3.16) 

The ie prescription is such that the time coordinate ti of operator O, is greater by —ie 
than the time coordinate of any operator to the right of Oi [Ml 135] . We compute the 
associated non-relativistic three-point correlator by Fourier-transforming with respect to 
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the Xj . This is done in appendix and one finds 

Gi23 - T{5,)T{5,)T{5,) ^""^^-''^ [^u-^^) K2-^e) 

x+ x + 



X 



2 ^13 iMoX^ /^23 f 



X (Ml - |M2|)^^+^^-ii?((5i, 52)$i(5i, 1 - (53, 5i + ^2, '^''7,^\ ^t;i2(Mi - IM2I)) 



I Mo 



(3.17) 



with 



^12 , ^23 -^13 



^12 = ^^^^i^ + ^^. P^. (3.18) 



^12 ^^ -^23 *^ -^13 *^ 



The functions S and $1 are defined in appendix \M This expression is consistent with 
(I3.10p as can be seen by noting that when x^ is positive, {x~^ — ie)~^ = |x'^|~''. The 
other Wightman three-point functions can of course be obtained by permuting the labels 
1,2,3. The relations given in (13. 9p can then simply be understood as particular examples 
of these permutations in which the labels 1 and 2 are swapped. These permutations are 
particularly simple because the coefficient in front of / in (13. 3p is symmetric with respect 
to the swapping of 1 and 2. Permutations of either 1 or 2 with 3 are less trivial because 
this coefficient does not exhibit symmetry under these swappings. These cases can be 
checked by using an integral expression for / given in eqn. (jA.Sp . 



3.4 Retarded and advanced three-point functions in momentum 

space 

For non-relativistic retarded and advanced three-point functions, we will apply the ap- 
proach used in section |2] for obtaining non-relativistic momentum space two-point func- 
tions directly from relativistic functions in momentum space. There, we saw that we can 
simply start with the relativistic function and rewrite the relativistic momenta in terms 
of non-relativistic momenta. The relativistic retarded three-point function was computed 
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in pgfl : 

^^'Vl,2;3) = ,.,.,. ' -, z^CT ^i)^Cr Piv)^^^ Pi) / dzz-^ 

X [-(wi - ze)2 + p{^^ + p2] "^1/2 K^^ [z^-{ujx - leY + p?^^ + ^2^ 

X [-(a;3 + lef + p^ + p^] ""''^ K.3 (^-Ic^s + ^e)^ + ^l^ + p^) . (3.19) 

The notation here is a straightforward generahzation of notation used in sections [2] and 
131 The cjj, "Pi^y and pi are relativistic energies and momenta. The spacetime is (i + 2 
dimensional, where one of the spatial coordinates, y, is singled out. The third operator 
is singled out in g\i{\-, 2; 3) because, in position space, this correlator corresponds to the 
case where the operator O3 has the largest time. The (5j and z/j are as defined in (13. 2 p and 
f l3.12p . Applying the map (12. lip to each of the Wj, 'Pi^y pairs gives 

GgVl,2;3) = ,,,,,, ,^ -, -bCr^,)bCy" Mi)b''Cy" V,) \ dzz2 

X [-4Mifii +p\ + ie{ni + Ml)] "'/' K,, (z^J-AM.ni + pj + ze(fii + Mi) j 

X [-4M2(]2 +pI + ^e(^2 + M2)] "'/' K,, rzy/-4M2fi2 + pI + ^£(^2 + ^'^2)] 

X [-4M3fi3 +pI- ieiQs + M3)] "'^^ K,, (z^-4M,n, + pj - ze{n, + M-s)) . 

(3.20) 

We will confirm in the next subsection that, in position space, G}j (1,2;3) is such that 
O3 has the largest time. The advanced non-relativistic three-point function is simply the 
complex conjugate of the retarded function: 

Gf = {Gfy. (3.21) 

These results will be used in the next subsection to compute retarded and advanced 
three-point correlators in position space. 



^In ref. }29j . the retarded function was computed in 3+1 dimensions in the case where all three 
scalar fields have the same conformal dimension, but it is easy to generalize this to arbitrary spacetime 
dimension and to three different conformal dimensions. 
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As in the case of the time-ordered three-point function, we can extract a bulk-to- 
boundary propagator: 

Qj^ ^ ^f +1 [-4Mn +p^- ieiyt + M)] "'^ K^ {z^~AMVl + p2 _ ie(yt + M)\ . (3.22) 

In inserting the factor za+i, we have again anticipated the AdS/CFT result. This is 
a retarded non-relativistic "bulk-to-boundary propagator". Notice that the three-point 
function G\^ contains two different types of propagators which are distinguished by the 
sign in front of the ie. We have identified the retarded bulk-to-boundary propagator by 
selecting the propagator which is analytic in the upper half il-plane. The other type of 
propagator appearing in G)j above is an advanced propagator: 

g^ ^ ^i+i [-4Mfi + p2 + ieiyt + M)] "'^ K^ (z^^AMnTf^^idSl^M)\ . (3.23) 

The fact that the retarded three-point function can be expressed as an integral over two 
advanced propagators and one retarded propagator was pointed out in the relativistic case 
in [22]. Here we see that this same structure carries over to the non-relativistic theories 
we are considering. In section HJ we will verify directly from non-relativistic AdS/CFT 
that it is appropriate to identify Qr and Qa as retarded and advanced bulk-to-boundary 
propagators by performing an analytic continuation from a Euclidean propagator. 
Note that the real-time propagators Qp-, Qr, and Qa satisfy the relation 

Gf = e{Q + M)gR + e{-Q - M)gA. (3.24) 

This relation is not surprising when we recall that the momentum-space two-point corre- 
lators satisfy a similar relation (I2.5ip . Eqn. (I3.24p can be rewritten in a form analogous 
to fl230D . 

gp = e{M)gR + ^(-m)^^, (3.25) 

as can be verified easily by rewriting the propagators as 

.2 \ ^/2 



gF ~ e~'^''^~''>'^\M['/^(-2n + ^-ieK'lS 



X K,[z^2\M\^~'^'^~''^'\I-'^^ + {^-^^m], (3.26) 



2 \ ^/2 



Qr = Q\ ~ e-'^-''>'^\M\''"(-2n + ^^-i^ 



X K, I z^/2\M\e'^'''-''^'\j -2n + ^-ie\. (3.27) 



23 



3.5 Retarded and advanced three-point functions in position 

space 

We will obtain the retarded and advanced three-point functions in position space by 
taking the Fourier transform of the momentum-space correlators computed in the previous 
subsection. As we saw, these correlators can be derived from relativistic three-point 
functions by redefining momentum space variables. The relativistic three-point function 
^]^ (1,2;3) which leads to the above expression for G*}^ (1,2;3) is the Fourier transform 
of the retarded three-point function for which the operator O3 has the largest time. It is 
natural to expect that the same will be true for the non-relativistic three-point function 
G}j (1,2;3), and this is indeed the case. For simplicity, we restrict ourselves to the case 
Ml > 0, M2 > 0, which in turn implies that M3 < 0. The details of the calculations are 
given in appendix [Dl For the retarded correlator, we find 

n(U \n(M ^r^^'>(^ 9-'^\ - (27r) t » 'SiJ^jMi) <+ \-S2\^+ \-Si\^+ \-53AMi(x^^+ie)/x+^ 

U[Mi)U[M2)Lrj^ [i,Z,6) ^ IX13I IX23I IX12I e ^13 13 

A f 

X B{6s, Si)<^ii6s, 1 - 62, 61 + ^3, -XT' -^^12^2) + ^(4i) 



X 



Mf^+'-^-^M^^-'B{63, S2)M^3, 1 - 5i, ^2 + ^3, -^, ivi2mA, 

(3.28) 



with 



2 2 

•^12 I -^23 "''IS 



^^^ ^ (3.29) 



"''12 "''23 "''13 



The function B{x,y) is the Euler beta function, while $1 is a hypergeometric function. 
Explicit definitions for these functions are given in appendix |Al Notice that the above 
expression for 9{Mi)9{M2)Gj^ is non-vanishing only when x^ > x^ and xjj" > xf, which 
is consistent with our expectation that G}j, (1, 2; 3) is a retarded correlator for which x'^ /2 
is the largest time. The advanced correlator, with x|^/2 now being the smallest time, is 
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l-'-lSl 1-^231 l"''12l 



r(5i)r(52)r(53) 



found to be 



9(M,)9{M,)Gf{l,2;3) 



X 5(53, 52)$i('53, 1 - 5i, 52 + ^3, -— i, z^;i2Mi) + ^(x+ ) 

IVI2 

X Mf -iM2'i+^^^-i5((53, (5i)$i((53, 1 - 62, 5i + ^3, -^, -^t;i2M2)|. 

(3.30) 

In position space, retarded and advanced three-point functions can be expressed in 
terms of the Wightman (non-time-ordered) functions [SB]. For instance, the retarded 
three-point function which corresponds to choosing x^ to be the largest time is given by 



.(3) 



G'^>{xt,xt;xt) = 9{xt,)9{xi,){Gs^2-G^s2 + G2^3-G 



231 y 



+ ^(42)^(4l)(G'321-G231+Gi23-Gi32). (3.31) 

We can use this formula to compute the retarded three-point correlator without making 
assumptions about the signs of the Mj. In the case Mi > 0, M2 > 0, it is easy to 
check using our formula for the Wightman function (I3.17P that only 6*123 and G213 are 
non- vanishing, and that 03.311) agrees precisely with 03.281) . 

We can find an alternate expression for the retarded three-point correlator in position 
space by taking the Fourier transform of the momentum space result 03. 201) : 



24-1 E. -«7r6-lr(Ai)r(A2)r(A 



G^(xt,X,M,^ - ' ^^' " ^U^lJil^2jil^3j,,>^,,, / ,. /,..,. ..-.-3 



^<5(^M,) f dz I dx^d'^xz- 



T{5,)T{52)T{5,)T{Y.A , 

X Qr{x - Xi,Z] Ml, /\i)<5r{x - X2, Z] M2, A2)Qr{x3 - X, Z] M3, A3), 

(3.32) 



where 



-z)^+i 



gR(x,z-M,A) = i-^|M|^-iz^^(a;+)|a;+|-^e'^^(-^+^^)/-^ (3.33) 

is the Fourier transform of the Schrodinger retarded bulk-to-boundary propagator 03.22p . 
This result assumes that Mi > 0, M2 > 0, M3 < 0. 
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4 Zero-temperature bulk-to-boundary propagators in 
momentum space 

We will now confirm that our proposals for the real-time bulk-to-boundary propagators, 
( ]3.13p . ( ]3.22p . and ( ]3.23p . are consistent with what are obtained from non-relativistic 
AdS/CFT. Since the results of the previous sections were derived without making recourse 
to AdS/CFT, it is reassuring that these results are reproduced when applying holographic 
duality prescriptions extended to non-relativistic theories [H [10]. In AdS/CFT, we may 
construct retarded propagators by first computing Euclidean propagators and then per- 
forming an analytic continuation [29]. To arrive at a Schrodinger- invariant theory, we 
must also perform (in the case of momentum space propagators or correlators) a change 
of momentum variables, as we saw in earlier sections. Once we have the retarded propa- 
gator, we may obtain the remaining real-time propagators from this by exploiting certain 
relations of which f l3.24p is an example. 

The zero-temperature Schrodinger metric in Euclidean space Poincare-like coordinates 
has the form [IH [HI [13] 



dsl = r\l + f3^r^)dT^ - 2il3\*dTdy + r\l - f3^r^)dy^ + r^dx^ + ^. (4.1) 



We are following the notation of Herzog, Rangamani, and Ross [TT]. The radial coordinate 
r ranges from zero to infinity, and the spacetime has a boundary at r = oo. The Euclidean 
time r and the Minkowski time t are related by r = —it. We lump together d of the spatial 
coordinates, Xi, ..., x^, into a c?- dimensional vector x with x-x = x^; as in previous sections, 
we reserve vector notation for this d-dimensional subspace. Although the metric (14. ip is 
only a solution to (a truncated version of) type IIB supergravity when d = 2, we will keep 
d general throughout this section to facilitate comparison with the results of previous 
sections. 

It is important to stress that t should not be identified with physical time in the 
Schrodinger field theory. The physical timq^ x"^/2 is obtained by switching to light cone 
coordinates: 

x^ = y±t. (4.2) 

The time t does admit the interpretation as the time coordinate of a parent theory that 
lives in one higher dimension relative to the Schrodinger theory we wish to study. As 
we described earlier, the parent theory is a non-local dipole theory [T71 [TH [19], and 
the Schrodinger theory is obtained from this by applying DLCQ to the x~ direction. 
Our approach will be to compute quantities in the parent theory and then perform a 



^This coordinate was denoted u in [TT]. 
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transformation (the momentum-space counterpart to fj4.2p ) to obtain quantities in the 
Schrodinger theory. 

We begin by computing the retarded propagator of a minimally coupled massive scalar 
(p in the parent theory. This can be obtained by first finding — in Euclidean signature — 
the solution ipE to the scalar wave equation which is well behaved in the interior of the 
spacetime. After Fourier-transforming along all directions except r, 

FE{ujE.PyJ,r) = j drdyd^xe-^'^^^-'Pyy-'^-^M'r.y.x.r), (4.3) 

the wave equation has the form: 

F'^ + ^Fi, -^A^I+pI+p" + (^"(Py + ^^^y^" + ^"^"] Fe = 0. (4.4) 

The prime denotes differentiation with respect to r, and m is the mass of the scalar. The 
solution to this equation which is well behaved as r — )■ is 

C 

FE{ujE,Py,p,r) = ^—K^{k/r), (4.5) 

with 



and 

u = y ^^^^ + m2 + /3^(py + iueY. (4.7) 

The constant C can be fixed by normalizing the solution at some large but finite value of 
r. If we now analytically continue ue -^ —i{uj + ie), then the resulting function will be 
analytic in the upper-half u plane and will thus admit the interpretation as a retarded 
propagator in the parent theory. 

We would like to re-interpret this result in terms of the Schrodinger theory. We do 
this by first inverting the coordinate transformation (14.21) : 

—i_ 1 _ 

r = — (x+-x"), y = -{x+ + x~). (4.8) 

Applying this coordinate transformation to (14. 3 p allows us to identify the frequency Q 
associated with the time coordinate x~^ and the momentum M associated with the coor- 
dinate x~: 

^ = ^{^-Py), M=-{uj + Py). (4.9) 
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The inverse of this map is 

uj = n + M, py = M-n. (4.10) 

These mappings are of course identical to (12.101) and (12. lip . The retarded propagator 
parameters can then be rewritten as 

k = ^/-AMn +p^-ie{n + M), (4.11) 



Setting r = 1/z, we then have 

Fe -^ Gr = Cz^+^K^ fz^/-AMn + p2 _ ie{n + M)) , (4.13) 

which is consistent with our earher finding for the retarded bulk-to-boundary propagator, 
eqn. (I3.22p . if we identify the conformal dimension A according to 



A = ^ + l + iy=^ + l + y'il±^ + m2 + 4/32M2. (4.14) 

We note the similarity between the AdS scalar bulk-to-boundary propagator and the 
scalar propagator given in (14. 5p . (14. 6p . and (14. 7p . with the main difference being the extra 
term /3^{py + iueY inside the square root in (14.70 . The AdS case corresponds, of course, 
to having /3 = 0. This change in v is the only effect of having a nonzero /3, at least as far 
as the scalar bulk-to-boundary propagators are concerned. 

We can apply this procedure to the other types of real-time bulk-to-boundary prop- 
agators identified in |29]. One finds for instance that the Feynman and advanced bulk- 
to-boundary propagators constructed in this way agree precisely with (IS.lSp and (13.231) . 
again writing the conformal dimension as in (I4.14p . It was argued in [29] that one can 
consistently define Wightman bulk-to-boundary propagators by extending the applicabil- 
ity of standard fiat space circling rules. This approach naturally commutes to the present 
context, with the Wightman propagators given by 

G+ = Gf-Qa. G' = Gf-Qr. (4.15) 

Combining these expressions with (I3.24p . it is straightforward to derive that the Wightman 
propagators can also be written as 

^^ = 2 [sgn(fi + M) ± 1] Im^R. (4.16) 
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It is easy to check that the various bulk-to-boundary propagators identified in this 
and the previous section reproduce the momentum-space two-point functions computed 
in section [2l This essentially follows from standard AdS/CFT two-point function com- 
putations. Since the momenta play a passive role in the AdS/CFT prescription which 
takes us from bulk-to-boundary propagator to two-point correlator, it makes no difference 
whether we apply the mapping f l4.10p before or after we execute this prescription. More- 
over, the different ie prescriptions which distinguish the various types of bulk-to-boundary 
propagator are also passive and carry over directly to the resulting two-point functions. 



5 Finite-temperature correlators 

5.1 Black hole metric and bulk-to-boundary propagators 

Finally, we turn to computing Schrodinger correlation functions at finite temperature from 
non-relativistic AdS/CFT. The previous sections have established the various ingredients 
we will need to obtain these correlators, in particular the momentum-space mapping 
(14.1 op from the parent to the Schrodinger theories and the identification of the different 
real-time propagators. We will begin by using these methods to construct the different 
real-time bulk-to-boundary propagators at finite temperature. We will first obtain the 
retarded propagator by computing a thermal Euclidean propagator and then analytically 
continuing this and using f l4.10p to obtain a thermal retarded propagator for the gravity 
dual of the Schrodinger theory. The remaining real-time propagators can be obtained 
from this using various identities borrowed from the case of a relativistic CFT. With 
these propagators in hand, we will then proceed to compute various two-point correlators 
as well as the retarded and time-ordered three-point correlators with the help of circling 
rules as in ^^. 

To study the theory at finite temperature, we consider the following (Euclidean) black 
hole geometry [HI [T21 [12] : 



dsl = r^k{r)-^^^ [-f3^r^f{r){idT + dyf + f{r)dr^ + dy^ + k{r)dx^] 
. dr"^ 



+ Krf"4^,. (5.1) 



% Hr)^l^^ 



/(r)^l-^, k{r)^l + ^. (5.2) 



r+ is the location of the event horizon. As in the case of zero temperature, we are adopting 
most of the notational conventions of [11]. In this section, we restrict attention to the case 
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d = 2 (i.e. X is two-dimensional). In order to simplify some of the expressions that will 
follow and to make contact with previous work, we will define a new radial coordinate u: 



«^^. (5.3) 



The metric then reads 



dsl = !±A:(w)-2/3r «2.2/H 

u 



(5^rl'-^^iidT + dyf + f{u)dT^ + dy"" + k{u)dx^ 
u 



/(n) = l-u2, k{u) = l + l3'^rlu. (5.5) 

The temperature can be computed from the surface gravity, which in turn is computed 
from the Killing vector generator of the event horizon, ^, according to the formula 

«' = -^(V'^e')(V.6). (5.6) 

The vector ^ is normalized so that the component along the Schrodinger time direction is 

1: 

d d d 

dt dx+ dx~ 

With this definition, we find that the temperature of the Schrodinger theory is given by 

twice that of the parent dipole theory, 

T5 = f = ^ = 2r„ (5.8) 

since Tp should be measured with respect to the time t. From (15. 7p . we see that we also 
have a chemical potential // associated with the conserved charge ^|^. With our set of 
conventions, 

/i = -2. (5.9) 



As was noted in |TTl[T2], these observations can also be discerned by considering a canon- 
ical ensemble for the parent theory with a density matrix given by 

where the Hamiltonian Hp = ^, Hs = S^pr, and P_ = g|^. We see that the thermal 
Schrodinger theory is naturally described in terms of a grand canonical ensemble, where 
P_ is interpreted as the particle number operator. 

30 



We will focus on the case of a minimally coupled massless scalar for simplicitjo. In 
momentum space, the equation of motion has the form: 



F'' 



u- 



u{l — u^) 



F'- 



2\'- E 



U^il — m' 



2^2 



ulu + {pi + p')u{l - V? 



where we have defined 

A = /3r+. 

We have also defined the dimensionless quantities uje-, Py, and p according to 

^E „ _ Py -_ P 



UiE 



Py 



P 



(5.11) 
(5.12) 

(5.13) 



2r+ ' ^ ' 2r+ ' ' 2r+ 
To solve the equation of motion, we use the following ansatz for the Euclidean propagator: 



FE{uJE,Py,P,u) 



U^(l-U 



i^sf^i 



l+Ur^'''H{u^E.Py.P 



U 



The resulting equation for H is the Heun equation. 



so long as we takq^ 



u u-1 u-d 



H' + 



u{u — l)(n — d) 



-H = 0, 



\ = l-Jl+k\Py 



iuje) 



(5.14) 



(5.15) 



(5.16) 



The Heun parameters are 



7 = l-2y'l + A2(pj, + iLL;E)2, 5 = 1 + a;^, e = 1 + icj^;, 
a = (3 = — 



uje + 1- Jl + k^{Py + ii^Ey, d 



q 



-UJe 



2Jl + A^{py + iu;Ey-l 



OJi 



pI-p'- 



(5.17) 



^ An example of such a scalar is the fluctuation in the metric component gi^'j^ , which has been shown 
to decouple from all other fluctuations [TTJ [T2] . As we show in appendix [El this scalar fluctuation can 
be made massive by considering non-trivial harmonics on the internal 5d space, in which case the wave 
equation generalizes to F'^ ~ J^S) Pe " „^(i-«^)^ | ^l"+ (P^ +p^ + |A^n^)M(l -u^) + [A'^{py + iu}E)^ + 

jl{l + 4)](1 — u^) >Fe = 0, where / and n^ are quantum numbers defined in the appendix. 

^°We have chosen the root which reproduces the zero-temperature limit, in which A — > and A ^- 0. 
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The Heun's function H must be chosen to be the solution which is well behaved near the 
origin n = 1 of the Euclidean spacetime, which for positive frequencies is o 

H{ujE,Py,P, u) = HI{1 — d, af3 — g; a, /3, 5, 7; 1 — u). (5.18) 

Analytically continuing cve -^ —i{Lo + ie), using fl4.10|]. and defining fl = r2/(2r+) and 
M = Af/(2r+), we find that the retarded propagator ia^n 

^R(ri,p,n;M) ~n\l-u)^5(«+A^)(l + M)i("+A^)i7(fi,p,M), (5.19) 

with 



A = 1 - V1 + 4M2A2, 7 = 1- 2V1 + 4M2A2, 6 = 1 - i{n + M + ie) , 

e = l + fl + M + ie, 

a = /3 = ^—^ (O + M + ie) + 1 - VlTiMV^, 

q = -i^ {n + M + le) h^lTAM^ ~ + ^^^ - p'^ + ie{n + M) . 

(5.20) 

Since it is not possible to normalize Qr aX u = 0, we introduce a cutoff mb ^ 1, and 
normalize the propagator at m = m^: 

As in [29], we define the thermal Feynman bulk-to-boundary propagator in the black 
hole background to be given by 

QF{^,Py,P,u) =ReQR{uj,Py,p,u) + icot\i{TTiv)lmgR{iv,Py,p,u). (5.22) 

Although [29] focused primarily on the gravity dual of thermal A/'=4 super Yang-Mills 
theory, the arguments (which were based on extending field theory finite-temperature 
circling rules to bulk-to-boundary propagators) were quite general and should apply to a 
wide class of theories including the dipole theory conjectured to be dual to the black hole 



^^See ref. |37| for a comprehensive discussion of the Heun equation and its solutions. A brief review of 
aspects essential for the present context is also given in appendix E of ^^. 

^^It is amusing to note that if we keep the ie in the exponent as well so that Gr ^ (^ — u)-5(^+-'^+"=)^ 
then the retarded propagator vanishes at the horizon m = 1 as opposed to the usual boundary condition 
that it behave like an in-coming wave. 
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spacetime fl5.4p . (Note that thermal A/'=4 SYM corresponds to the special case A = 0.) 
Applying f l4.10p . we obtain a similar relation for the Schrodinger propagators as well: 

gpi^, p, u; M) = ReGRin, p, u; M) + i coth (7r(rj + M)) Im^i? (rj, p, u; M). (5.23) 

This formula allows us to compute the Feynman bulk-to-boundary propagator directly 
from the retarded propagator we have already calculated. Again borrowing results from 
the relativistic case, we may compute Wightman propagators using the same relations we 
have at zero temperature |29j : 

G+ = Gf- Ga, G~ = Gf- Gr. (5.24) 

In conjuction with (I5.23p . these expressions imply that 

\m.GR. (5.25) 



X _ gT27r(f2+M) 



5.2 Real-time thermal two-point correlators 

The real-time bulk-to-boundary propagators obtained in the previous subsection can be 
used to compute real-time two-point correlators at finite temperature. We start with the 
retarded two-point function, which arises from a boundary term in the quadratic part of 
the effective 5d actioio of the form J dD.(Ppy/^gg'^^(l){—i^,—p)du(p{i^,p}\u=UB- Setting 
the scalar field (j){Q,p, u) = Gr{^,P, u)(J)b{^i'P) and varying with respect to the boundary 
value (pB, we find 

Gr{SI,p-M) = -2Ny^g'^'^VR{n,p,u-M)duGR{n,p,u-M)l^^^. (5.26) 

iV is a supergravity normalization factor. In order to evaluate this expression, we must 
expand the Heun's function H near the boundary ai u = ub- This is facilitated by first 
expressing the Heun's function (I5.18P as a linear combination of independent solutions 
defined in the vicinity of u = 0: 

H{u) = H^{u) + AH2{u), (5.27) 

where 

Hi{u) = Hl{d,q,a,P,-f,5;u), (5.28) 



^•^ We may take the effective 5d action to be any of a number of consistent truncations of IIB supergravity 
that have been found in the hterature (TTJ [T^J [131 [3H1 [Ml [101 [H] • Since we are considering a generic 
minimahy coupled scalar, our results are independent of this choice. 
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and 

H2{u) = u^-''Hl{d, g - (7 - 1)(5J + e), /3 - 7 + 1, a - 7 + 1, 2 - 7, 5; u). (5.29) 

The coefficient A is determined by matching flS.lSp and fl5.27p at some value of u between 
and 1; in practice, this must be done numerically. Since our formula for the two-point 
function, eqn. f l5.26p . only depends on the behavior of Qr and its ffist derivative near 
ti = Mb <S 1, it suffices to expand Hi and H2 about u = 0: 

Hi{u) = l + ciu + C2U^ + ..., H2 = u^ + ..., (5.30) 

q q{e - S + q) + a^-f 
ci = , C2 = -— — — T , 5.31 

7 27(7 + 1) 

where we have set the Heun parameter d = —1. After throwing away terms which diverge 
in the n^ — )■ limit, we find 



witfJi^l 



GR{n,p-M)=4riN 



X-2{A + C2) + cl + -^{n + M) 



(5.32) 



We can obtain an analytic expression for A in the limit of small fl, M, and p by plugging 
the ansatz H{u) = l+qh{u) into (I5.15P and expanding the Heun parameters. The solution 
which is regular at m = 1, namely h{u) = log(l + u), gives a solution H{u) which is valid 
for all u so long as we only keep the leading order terms in q. Expanding this solution to 
second order in u about u = then gives 

A+C2 = ^^(^ + M) + ip2 + o{fi^, fiM, M^ op^ Mp^p^). (5.33) 

Plugging this into (I5.32p yields the thermal retarded two-point correlator in the limit of 
small frequency and mass. This quantity was also obtained in [12] for the case p = 0. 

To compute the Feynman two-point function, we can perform a similar analysis em- 
ploying the various thermal bulk-to-boundary propagators along with finite temperature 
circling rules. We will now give a brief review of circling rules as they apply to Witten 
diagrams. Each vertex of the diagram can be either circled or uncircled, and a general 
real-time correlator will be given by a sum of diagrams having particular vertices circled. 
The circlings determine the type of real-time correlator (e.g. time-ordered, retarded, etc.): 



^^Note that if 7 is a negative integer, then Hi can no longer be expanded in a simple power series — 
logarithmic factors will also appear. See appendix E of [29j for a more complete discussion of this 
subtlety. 
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if Qi extends between two uncircled vertices, then Qi = Qp-, the Feynman propagator. In 
the case of two circled vertices, we instead have Qi = Qp = —Q*f- Since the momenta pi 
are defined to all be out-going, Qi = Q~ when the bulk vertex is circled and the boundary 
vertex is uncircled, and Qi = Q^ when the reverse is true. If the number of circled ver- 
tices in a given diagram is odd, then the contribution from that diagram also receives an 
overall minus sign. In particular, for the time-ordered two-point function, the vertices at 
the boundary are uncircled, and the derivative-coupling bulk vertex can be either circled 
or un-circled. Summing over both contributions yields 



GF{n,p;M) = -2Ny^g''''iQp{n,p,u;M)duQFi^,P,u;M) 

- Q~{n,p,u;M)duQ-^{n,p,u;M)] , (5.34) 

/ u=ub 

which, upon substituting fl5.23p and fl5.25p . leads to 

Gf{'^, P; M) = ReGR{n, p;M) + i coth (7r(ri + M)) IraGnin, p; M). (5.35) 

Recalling that fi = -^ and M = ■^, where fi and M are the physical (dimensionful) 
frequency and mass, we see that this is precisely the standard Kallen-Lehmann relation at 
finite temperature with energy 2Q, particle number M, and chemical potential fi = —2. It 
is straightforward to check that in the zero-temperature limit where r_|_ — )■ 0, this reduces 
to the relation we found in section [2^ eqn. (I2.5ip . 

A similar application of circling rules in the case of the Wightman functions leads to 

5.3 Thermal retarded and advanced three-point correlators 

The retarded three-point scalar correlator (which is equal to the sum of all diagrams with 
vertices either circled or uncircled with the exception of the vertex with the largest time 
which remains uncircled) has the same universal expression in the black hole background 
which is dual to the dipole theory as the one obtained in 



9r\pi,P2;P3) ^ 5'^iPl+P2+P3) / du^/^Q*J^{pl,u)Q*Ji{p2,u)QR{ps,u). (5.37) 

JuB 

The position-space version of this correlator is such that the operator O3 has the largest 
time. Note that the four-momenta pi are defined to be all out-going, i.e. they point toward 
the boundary. This result can be obtained in a number of ways, including the application 
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of circling rules to Witten diagrams (as advertised) or by analytically continuing the 
Euclidean three-point function. 

At this point, we should make a few comments regarding the integration region in 
fl5.37p . In the case of thermal A/'=4 SYM, it has been proposed that one must consider 
the maximally extended AdS-Schwarzschild black hole spacetime in order to reproduce 
Schwinger-Keldysh correlators from AdS/CFT [2B]. That is, in addition to the quadrant 
which is integrated over in fl5.37p . there are three more quadrants that one should also 
take into consideration. However, as is argued in [29], the use of circling rules renders 
unnecessary the introduction of additional quadrants. This result is useful since it avoids 
one of the more subtle issues which arise in real-time AdS/CFT calculations, namely the 
role of global properties of the bulk spacetime. We will keep exploiting this simplification. 
It was also pointed out in [25] that in the case of a retarded three-point correlator, 
neglecting extra quadrants is consistent with the expectation that one should need to 
integrate only over a maximally causal region of the spacetime that includes the original 
boundary at u=ub- (In [29], mb=0.) This is tantamount to integrating over only the 
original quadrant — from the event horizon at u=l to the boundary at u=ub — since the 
other three are causally disconnected from it, as is readily apparent from the fact that 
the remaining quadrants lie beyond the event horizon of the black hole. ( For a different 
perspective reaching the same conclusion see |42j.) 

Having argued for the validity of the three-point function (15.371) for the parent theory, 
it is then a simple matter to transform this to a Schrodinger three-point function using 
(I4.10p for each of the three momenta: 

G%\n,,p,-M,) ~ 5(^M,)5(5^a,)5(^p,) 

i i i 

X I du./^gu^,,pi,u-M^)g*^{n2,P2,u;M2)gR{n3,P3,u-M3). 

Jub 

(5.38) 

The thermal retarded bulk-to-boundary propagator Q^ was given in (15.211) . The thermal 
advanced three-point function is given by the complex conjugate of (15.380 . 



5.4 Thermal time-ordered three-point correlator 

In this section, we compute the Schrodinger thermal time-ordered three-point correlator. 
As we have done in the case of the retarded three-point in the previous subsection, we 
could first compute the time-ordered three-point function in the parent theory and then 
transform this using (I4.10p to obtain the Schrodinger time-ordered function. However, 
considering the various results we have obtained in earlier sections, it is clear that this will 
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p,^ ' 




p,^ 



u=l u=0 u=l u=0 

Figure 1: Witten diagrams contributing to the time-ordered three-point correlator. 



be equivalent to a direct computation using a non-relativistic version of Witten diagrams 
and circling rules. A generic three-point Witten diagram is of the form 



«s 



(5.39) 
where the types of bulk-to-boundary propagators Qi appearing in fl5.39p depend on which 
vertices are circled. 

In the case of a time-ordered function, we must add the contributions of all diagrams 
which result from considering all possible circlings of internal vertices, with the external 
vertices uncircled. For a three-point function tree-level diagram there is only one inter- 
nal vertex, and the two contributing diagrams are shown in Fig. 1. A straightforward 
application of the circling rules reviewed above gives 



GP{n,,p,) ~ 6{J2M,)d{J2^i)6{Y, 



Pi 



du^/^QF{^l,Pl,U]Mi)QF{^2,P2,U]M2)QF{^i,PZ,U]M2, 



ub 

1 



du^/^g {Qi,pi,u;Mi)g {Q2,P2,u]M2)g {Vt-i,p^,u]M^ 



Ub 



(5.40) 



This result is consistent with (13. lip in the zero-temperature limit since in this limit, the 
second term vanishes. This can be seen from (I5.25P by writing Q, + M = ^^^ and sending 
r+ — !■ while keeping Q and M fixed. Momentum conservation ensures that Qi + Mi < 
for at least one of the three sets of momenta, and the corresponding Wightman propagator 
Q~{Qi,Pi,u; Mi) will therefore vanish. 
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A Zero- temperature time-ordered three-point func- 
tions in position space 

We begin with the relativistic conformal three-point correlator: 

gP = -{0\TO^ix^)O2ix2)O,{x,)\0) = 
1 1 



F23^23 + ^23 + ^^J ^ 



(A.l) 



[^r3^i3 + 2^13 + ^^r i^uXu + ^2 + ^^r 

where Xij = Xi — Xj. Fourier-transforming with respect to x^ and introducing Schwinger 
parameters gives 

gP = --}—l——— / ds,ds2dsss'r's'rht'-' / dx^dx^dx^e-'^'^-^-^''^''^-'''^^^ 
T{di)T{d2)T{d3) Jo J 

It is useful to change integration variables from x^ to Xi^, X23, and x^ . The x^ integration 
then gives a (5-function of the form 6{Mi+M2+M3). The Xi^ and 0:^3 integrals also produce 
(5-functions, and we obtain 



r(50r(52)r(53 



GP = - '^W.r^...Z)'r' '^ ds,ds2ds,st~'sr'4-' 



3 

2 



(A.3) 

Evaluating the 5-functions, we find 






-|X'23| 'p'lal ~n '" "■ " """e ~" '" " "" 

(A.4) 
where 

/ ^ rd.3^rW2+X^2^3|^-Wl-^^2^3|^-^e--^^-^ ( M,-XX,S, \ ^ /M^+^\ 
Jo \ ^13 / \ ^23 / 

(A.5) 
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and 



^12 = + — + + + — • (A.6 

•^12 "''23 ""-IS 

In order to perform the integral /, it helps to consider separately the eight different cases 
which arise when we suppose that each of x^2i ^135 ^^^ x^g have definite sign. We adopt 
the notation /jj^, where i, j, /c = ±, such that lijk = I for the particular choices of sign 
specified by the i,j, k. For example /++_ is equal to / when xf2 > 0) ^13 > ^^^ x^g < 0. 
We need not consider /+_+ and /-+- since the corresponding inequalities are not self- 
consistent in these cases. Furthermore, a careful inspection of flA.SP reveals that it suffices 
to compute only the three integrals I+jk since the I-jk can be expressed in terms of these. 
In particular, we have 

/_++(5i, 52,^3, Ml, M2,W12) = /+++(52,5i,53,M2,Mi,-t;i2), 
/__+(5i,52,(53,Mi,M2,Wl2) = /++_(52, (51,^3,^2, Ml, -t;i2), 

/___(5i,52,53,Mi,M2,t;i2) = /+__(<52,(5i,53,M2,Mi,-i;i2). (A.7) 

In terms of these quantities, I is given by 

I = e{xtMxt3nxt3)h++ + ^(X+ )e(x+3)e(x+ )/++_ + ^(X+ )^(4i)e(42)/+„_ 

+ ^(4i)^(43)^(43)^-++ + e{xt,)9{4inxt3)i-+ + ^(x+i)e(4i)^(x+ )/___. 

(A.8) 
We begin by computing /+++: 

POO 

I+++ = \xt2\~^' / dzz^''-\M2 + ^)^^-^(Mi - zY^-^e'^^'^'eiMi - z)e{M2 + z). (A.9) 



We arrived at this expression by introducing a new integration variable via S3 = x^2-^- To 
proceed further, it is useful to note that the step functions impose the constraints Mi > 
and Ml + M2 > 0. We may then break up the integral as follows: 

h++ = |x+ |-^«^(Mi)[0(M2)/S+ + ^(-M2)^(Mi - |M2|)/S+]. (A.IO) 

We first focus on /+++. When M2 > 0, the z integration limits are and Mi. Taking 
z — )■ Miz, we find 

J«^ = Ml^+^^-^ M^2"'^ I dzZ^-'-\l - zf-'^^ (l + ^^ ' g-i2Mi._ (^_^^^) 

We can evaluate this integral by making use of the identity |13] 

t/xx''-^(l-x)^~^(l-/3x)-^e^^ = 5(i^,A)<I)i(z/,p,A + i/,/3,/i), (A.12) 
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where B is the Euler beta function, 



B{x,y)= f dtt^-\l-t)y-\ (A.13) 

Jo 



and $1 is a hypergeometric function of two variables 

(Q)m.-I 



$i(«, /3, 7, X, y) = f; [^\-^-^^K -y-. (A.14) 



m,n=0 



This identity holds provided Re(A) > 0, Re(i/) > 0, and |arg(l — /3)| < tt. For our 
purposes, the first two conditions will be always be satisfied if 6i > 0, while the second 
condition requires /3 = —M1/M2 < 1, which is true in the case of the integral /+++, for 
which both Mi and M2 are positive. We find 

/S+ = Mt+'-'-'M^'-'B{63, S2)<l>i{Ss, 1 - 61, 62 + ^3, -^,^Vl2Ml). (A.15) 

ik/2 

(2) 

The second piece of the integral, /+++, has the form 

/S+ = Mf2+^^-^|M2|^^-W dzz^'-\l-zf'-U-l + —^z] e™i^*^i^ (A.16) 

JlJM. \ |iw2| / 



We can again evaluate this by making use of the identity flA.12p if we first perform the 
change of variable ^ = (1 - \M2\/Mi)z + jMsl/Mi. We obtain 



/;{+ = {Mi-\M2\Y'+^^~^\M2f'~h''"^^^'^^^B{6i,6: 



2) 



X $1(^1, 1-^3.^1 + ^2, '^! '|^\ ^^12(Mi-|M2|)). (A.17) 

\Mo\ 



Putting the results for /i++ and /il+ together gives 



/+++ = \xt2\-''e{MAe{M2)Mi^^'^'-^Mi^-'B{5:,, 52)$i((53, 1 - 5i, 82 + ^3, -^, zt^i^Mi) 
+ e{-M2)9{Mi - |M2|)(Mi - |M2|)'^^+'^^-^|M2|^^-^e™i2l^^2l 
X 5((5i,52)<l'i(5i,l-53,5i + 52,^^%^/^,^t;i2(Mi-|M2|))j. (A.18) 

\M2\ ) 
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Similar methods lead to 

/"OO 

/++_ = Ix^sT"^' / dzz^-'-^\M2 + z\^^-^\Mi - z\^^-^e'''^^'d{Mi - z)e{-M2 - z) 
Jo 



X <l>i((53, 1 - 62, 61 + 53, ^, ivu\M2\) + e{\M2\ - Mi)M^'+''-'\M2\''~' 

IVl-l 

X 5(53,<52)$l(<53,l-5i,<52 + 53,^,^t^l2Mi)l. (A.19) 



/■oo 

^•■^ / dzz^^-^\M2 - z\^'-^\Mi + 2|''2-^e-*'^^2^^(Mi + z)e{M2 - z) 
Jo 



fOO 

1-53 / 
'12 



kl2r'^'^(-^2)i ^(Mi)^(|M2| - Mi)e™i2^iMf^-\|M2| - M^'^^+'^^-i 



X 5(^2, <5i)$i(<52, 1 - 53, 6i + 62, ^' J^'[ ^v,2i\M2\ - M,)) 

+ e{-M^)\M,\'^-^\M2\''^'■'^'B{5:,,5^)^,{5,^ - 82^5, + 5-,,-^^,tv,2\M2\)\. 

(A.20) 

B Zero-temperature time-ordered three-point func- 
tions in momentum space 

To compute the time-ordered three-point function in momentum space, we will Fourier- 
transform the position space correlator. It turns out that the most efficient way to perform 
the Fourier transform is to start with the relativistic function in position space expressed 
in terms of light cone coordinates: 

gf = -(0|rOi(a:i)O2(a;2)O3(x3)|0)= ^ 



F23^23 + ^23 + '^^\ ^ 
1 1 



[x^3Xi3 + xfg + ieY^ [xt2Xi2 + XI2 + «e]^3 



(B.l) 



We could start directly with the non-relativistic result f l3.3p . but this approach proves 
more cumbersome. Taking the Fourier transform with respect to all coordinates and 
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introducing Schwinger parameters leads to 

\d+2(_-\5i+52+S3 



5i-l„52-l„<53-l 
3 



~^o^ (0Tr\d+2( ■\5i +62+63 foo 

The momentum variables f2j, M, and pi are conjugate to x^, x^ and Xj respectively. To 
obtain this expression, we have switched integration variables from Xi to X13, X23 and X3, 
and we subsequently computed the d + 2 integrals over X3, giving rise to the ^-function 
factors out front. The integrals over xj"3 and x^3 give rise to 5-functions, and the integrals 
over X13 and X23 can be evaluated using 

^d id i{s2+S3)x\^+i{s2+S3)x1^-2is3Xi3-X23-ipi-xi3-iP2-^23 

(^TT\d i °lPl+°2P2+^3P3 

e * »l»2 + si«3 + ''2''3 . (B-3) 



(S1S2 + S1S3 + S2S3)'^/2 

We then find 

j-,-1 j^-i ^s^,-i , np!+«2Pi+^3Pi /■ ^ ^ -o +^0 + 

X -777:6 ■* »i»2+=i»3+''2''3 / (iXiqdXgqe ^""is+^'axaa 

(S1S2 + S1S3 + 5253)^^/2 y 13 23 

X (5((s2 + S3)x+3-S3X+3-Mi)5(-S3X+3+(si + S3)x+3-M2). (B.4) 

Using the (5-functions to evaluate the integrals over x^.^ and x^3 yields 

S^^-'-S^^^-'-So^""^ , (4A/ini-pf)si + (4M2n2-P2)^2 + (4M3n3-p§)s3 

X =^ = 2 6* ''1^2 + ''l''3 + ''2^3 . (B.S) 

(S1S2+ S1S3+ 5253)2+^ 

To facilitate evaluation of the Sj integrals, we make the following change of variables: 

S1S2 + S1S3 + S2S3 
«i = , (B.6) 
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under which the integral becomes 

/■oo 

X / n^"i<'"'^'^'*'''''"'''''""^^'''("i+"2 + W3)"^'''+^+'. (B.7) 

Jo ^. 

We have defined 

d ^ . d 



^,=^5.---l = A,---l. (B.8) 

Further progress is made with the help of an additional Schwinger parameter: 



/•CXD /"OO 

X / dzz^^^'^-''-' J]dn,n7-'e3[4^'^^^^-P?+-]/-^e-^'"^ (B.9) 



The integral over each Uj can now be identified as being essentially a modified Bessel 
function: 

gP = ,,,,,, ,^ — , — -6(yQ,)6(yM,)6''(yp,) 

X /" dzz-^ JJ(-4M,^- + p2 _ leY^/'^K,^ (^z^ -AM^Vt^ + p] - it^ . (B.IO) 
Jo J. 

In this expression, we have rotated the 2;-integration contour to eliminate a factor of yi 
that would otherwise appear in the argument of the Bessel function. This rotation can be 
done safely since the branch cut of the Bessel function lies along the negative real axis, 
and the square root in its argument is defined to be such that its real part is non-negative. 
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C Zero-temperature Wightman three-point functions 
in position space 

We begin with one of the relativistic Wightman functions in position space: 

1 
5'i23(2:i,a;2,X3) = -(0|Oi(xi)O2(a;2)O3(x3)|0) = 



F23^23 + ^23 + ^^(,^23 ^23)J ^ 
1 



X 



Fl3^13 + ^13 + ^^(^13 ^13)] ^ Fl2^12 + ^12 + ^^(^12 ^12)1 ^ 



(C.l) 



We factor out Xij—ie from each of the three denominators and take the Fourier transform: 



G 



123 



-5i/^+ 



Xno — 26 MX10— «e "xVn— «e 



-23 



'13 



S2 I 



12 



-<53 



f dx^dxl^dx^e-'^'^""^ "^^'^^^^ -^m^x^ 



X 



X 



^23 + 



2 + 
"''23''' 23 



[X 



\2 I ^2 



X 



12 



We have defined 



ei =e 



23 yl 

2 + 
•'''12"''12 

(X+ )2 + 62 



^T+ "^2 I ™2 
^"''23/' "T "''23 



n 


-<5i 


261 




- 


-<53 


^£3 





2 ^+ 



Jy 1 Q iL 



X 



13"^13 



13 



X 



13; 



\2 I ^2 



«e2 



-S2 



(C.2) 



X 



23/1 



\2 I ^2 



62 = 6 



^T.+ ^^2 I ™2 
^"''ISJ "T "''13 



[X 



13; 



\2 I ^2 



63 = 6- 



("''12/ ~r "''12 



[X 



12; 



\2 I ^2 



(C.3) 



Since the 6j are positive-definite, we may introduce Schwinger parameters. After changing 
integration variables from x~ to x~[^, X23 and x^ , we arrive at 



G 



27r{-iY^+^^+^^{J2iM, 



123 



i ^> ( ^+ 



V{5^)V{52)V{5^) 



(x 



23 



ieV^Hxt - ler^Hxto - i^ 



'13 



12 



dsids2ds:isl'-\l^-\sl'-^ I dxr.rfx7.e-^*^^^i3-*A^2X2-3 



'13""'' 23"= 



-si 



X e 



fi _,>- ,- ^23^23 
"1 *^23 \,+ )2+,2 



-S2 



^^— 13-^^^^^y^ 



-S3 



e\—ix^r)—i 



12 "(42)'+^' 



(C.4) 



The integrations on x^3 and x^3 yield (5-functions of the form 5{—Mi + S2 + S3) and 
(5(— M2 + Si — S3). Since the Sj only assume positive values, the first of these 5-functions 
imposes the condition that 6*123 will only be non- vanishing if Mi > 0. This is in keeping 
with expectations from perturbation theory, where a negative Mi would imply that Oi(xi) 
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is comprised of creation operators. Since the creation operators will annihilate the bra 
vacuum state (0|, G123 must vanish in this case. Furthermore, the two 5- functions together 
imply that the Sj integrals are supported on the curve si + S2 = Mi + M2. The factor 
S{Mi + M2 + M3) appearing in G123 then ensures that M3 = —Mi — M2 < when G123 is 
non- vanishing. This is once again in keeping with perturbation theory if we now consider 
Oslxs) acting on the ket |0). We also note that Mi + M2 > implies that when M2 < 0, 
we must have Mi > IM2I. These observations are useful for simplifying the remainder of 
the calculation of the Wightman function G123, as we will see. 

Evaluating the (5-functions arising from the ^^3 and X23 integrations gives 



r{Si)r{62)T{63 



<^123 t^/j: \-n/x \T-/i: \ 1^23 ^ *^J K^l?,-^^) {^12--^^) 



X e '''w' +' e ^"^23' +^ I[vi2)., (CS) 



where we have defined 



/■oo 

/>i2)= / ds,s'i-\M2 + s,Y^-\Mi-s-it-^e'^^^''e{Mi-s,)e{M2 + S3), (C.6) 



and 



2 + 2 4- 2 + 

_ "''12"''12 I "''23"''23 "''13"''13 , -(^ ^ \ ^ \ 

rft^ rv>^ ry^ 

^ "''12 I "''23 "''13 



- + ^^. fi^. (C.7) 



X12 ^e 3^23 ^^ "''IS *^ 



We have already evaluated the integral /(W12) in appendix |X1 It is given by 

A^i2) = \xX2\'n^^^ = Q{hh)\Q{M2)M{^^'^-^Ml^~^B{h,b2) 

X $1(^3, 1 - h. S2 + ^3, -^, tvuMi) + e{-M2)e{Mi - |M2|)e™-l^^^l|M2|^^-i 

X (Ml - |M2|)'^^+'^^-^E((5i,^2)$l(^l,l-(53,^l+(52 j^''~,^\ 2^12(Mi - |M2|))|. 

l-lVhi J 

(C.8) 
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Putting everything together, we find 

Gi23 Ti6,)Ti6,n6s) ^''''~''' ^""'''"'^ ^''''~''' 



X $i(<53, 1 - 5i, 52 + 53, -^, zz;i2Mi) + e(-M2)^(Mi - |M2|)e"'-l^^^l|M2|'^^-i 

M2 

X (Mi-|M2|)^^+^^^^i?((5i,52)$l(5i,l-53,5i + 52, '^''~,^\ 2t;i2(Mi-|M2|))j. 

IM2I J 

(C.9) 

D Zero-temperature retarded and advanced three- 
point functions in position space 

We will compute the retarded three-point function by taking the Fourier transform of the 
momentum space result we obtained in section 13.41 The retarded three-point function in 
momentum space can be written in the form 



r(5i)r(52)r(53)r(E.5,-f-i) ^^ '^ ^^ '' ^V -^o 

(2 \u\/2/ 2 \ ' 

2Mi y V 2M2 ) 



,2 \ ^s/2 

, i 
'3 



2M3 



X ( -2^3 + ^3 - ie I ^~i^e{~Ml)vi/2^~iTxe{~M2)v2/2^i-ne(~Kh)v3/2 



Pi 



X K., I z^2\M;\e-^-'^-''^y\\-2n, + ^ + ^^ 



P2 



X K., I ^V'2^e-^(-^^)/^/-2fi2 + ^ + ^^ 



pi 



X K.3 I ^v/2Me^^'^-''^^/^/-2f^3 + ^ - ^H • (D-l) 
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For simplicity, we will focus on the case Mi, M2 > and M3 < 0. Rotating the contour 
of the 2;-integration such that z — )■ 2;a/— z, writing the Bessel functions in terms of the 
following integral representation: 



V^ X*^ / , r / I I / "P^ 



/"OO i\M\ I „n I P^ f • 1 2 

Jo 
and performing the integral over z leads to 

B{M,)e(M,)Gf = r(^or(^2)rfe) ^(I:".)^'(I:''^.)^'(I:p") 

/■oo 

^ ^^(-2^.+4r+-)^^(-2n.+J,+.)^^ (-2^3+4-.) _ ^^_g^ 
We want to compute the Fourier transform of this expressionist 

Gf = ^^^^ j n dQd^e-^^'^^t+n-^.Gf, (D.4) 

One of the fij integrals (^3 say) eliminates the factor of 5(Xlj^i) appearing in Gp^. The 
remaining two f2j integrals then produce (5-functions. Similarly, one of the pi integrations 
removes the 5(^^pj), and the other two pi integrals are easily computed, with the result 

^ > \^-^l ^) I J„, J„, J„, „ o.\ „., ^ ,/„ ,/„ 



2 V "1 t»2 "3 . 



{u^U2Y'^-\u^u^y--\u2U;Y^'\ [ M^ M^ +\ / M2 M3 + 

X - , Y^ r _i I r X13 I I H 3^23 

(til + M2 + Ms)^" ' V Ml M3 J \ U2 Us 



X 



Further progress is facilitated by the following change of integration variables: 

1 U1U2U3 S1S2 + S1S3 + S2S3 
Si = ■ ■ , Ui = . (D.6) 

Ui U1 + U2 + M3 Si 



^^We include a factor of 2 for each of the three fli integrals because the non-relativistic energies are 
given by 217^. 
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After using the two 5- functions to perform the integrations over si and S2, we are left 
witfS: 

/■CXD 

where 



^23 / V ^13 



"''12 , "''23 "''13 



^12 = ^ + ^ - ^. (D.8) 



"''12 "''23 "''13 



Comparing the integral in (1D.7P with ( ]A.5p . we see that we have essentially already 



evaluated this integral. It can be expressed in terms of the functions /+++, I++-1 and J+ 

we defined in section ISTTj and appendix lAl We must first send Mi — )• —Mi, M2 — )■ — M2, 
and W12 — )■ — t'i2 and then impose the constraints Mi > 0, M2 > 0. From (IA.18I) . (IA.19I) . 

and (IA.20I) . we see that only /+ is non- vanishing under these constraints. Since / 

can be obtained from /^ by swapping parameters (see (1A.7I) ). this function will also 

contribute, and we find 

U{M,)U{M^)L. ^ - ^^^^^^^^^^^^^^^ |x,3| IX23I \x,,\ e 

X e^*^2(^23+i^)/43^(Mi)^(M2)^(4i)^(a;+)|^(x^2)^A^'"^^l'^'^'"^ 



M2 

X 5(^3, (5i)$i((53, 1 - ^2, 5i + (53, -— -, -W12M2) 

+ ^(X2+ )Mf^+'^^-lM2^^-li?((53, 52)<l>l(53, 1 - 5i, ^2 + ^3, -^, Zi;i2Mi) |. 

(D.9) 

Notice that this expression is non- vanishing only when xt, > X2 and x^ > xf, which is 
consistent with our expectation that G\^ is a retarded correlator for which x^ /2 is the 
largest time. 



16 



The e in (jD.5[) can be set to zero since the integral remains convergent in this hmit. A different e 



regulator is introduced in (|D.7p to impose the vanishing of G]^ at x^^ ~ and x^g = 0. The fact that 



G]j vanishes at these points is evident from (JD.SP since one or both of the (5- functions vanish identically. 
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The advanced three-point function can be obtained by starting from G\ = (G)^)*. 
We again suppose that Mi,M2 > 0, M3 < 0. Similar steps to the above lead to the 
expression 



12«3 



This time, the integral over S3 is precisely the same as / as it is given in flA.Sp . The 
conditions Mi > and M2 > now single out /+++ and /_++ as the sole contributors, 
and we find 



'^ - r(50r(52)r(53) l"^^^' ^^''^ ^""''^ 



X e^^2(^23+^^)/43^(Mi)^(M2)^(a;+3)^(x+3)<^ 0(x+)Mf^+''='-^M|^-^ 



Ml 

X B{63, (52)$l(53, 1 - Si, 62 + ^3, -T^T' ^^12^1) 

ikZ2 

+ e{xti)M('-'M^^+'-'-'B{5s, 5i)<^i{5s, 1 - ^2, 5^ + ^3, -^, -zt;i2M2; 

(D.ll^ 

We see that this result is consistent with the interpretation as an advanced correlator, 
with xf/2 now being the smallest time. 



E Sgi"^ metric fluctuation spectrum 

Let us begin with a brief reminder of the scalar spherical harmonics on 5*^, written as an 
S^ fibration over CP^: 

(is^s = (dx + sin^ /i co^Y + [dfx^ + sin^ l^i^i + 1^2 + cos^ /i cul)] 

Ui = - {cos ip d6 + sin ip sin 6 d(f)) 

a;2 = = - (sin ip dO — cos ip sin 6 dcp) 

CJ3 = -{dip + cos 9 d(l)) (E.l) 
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where < x < Svr, < ^u < f , < ^^ < 47r, < ^ < tt and < < 27r. 
The scalar spherical harmonics are solutions to the box equation 



Ds^Y 



1 
V9 



where 
and 



d.^g'^djY = -l{l + 4)F, / = 0, 1, 2 , 



Y = e*''^'^e^"<^'^e^"'*^e(^)C/(/x) 



(E.2) 
(E.3) 



cos I 



il"*-^0l(l + COS e)^>^+''^^2Fi ( - ne, ng + 2\n^\ + 1; 1 + |n^ + n^\ 



1 + cos 6^ 



U{fi) = (sin/i) 



2|n^|+2ne, 



COS/i) 



\n^—2n,i 



2-Fi(-n^,/ -nf, + 2;l + \n^- 2ra^|;cos^/i) 



/ 



n^ 



2n^\ + 2\n^\ + 2ne + 2n^ 



and where n^ < \n^\. 

In the black hole background p!T| [T2| [T3] 






A;(r) 



2/ /3'rV(r),^, , ^ ,2 fir), ,2 , rfy' 



r I — 

2 ,. , . ^2 



Mr) <* + "«> -iM^n-M 



dx' 



dr^ 



r'^fir) 



+ 



{dx + sin ULUsY 
k{r) 



+ ds. 



CP'^ 



(E.5) 



one of the simplest fluctuations, which decouples from the rest, is the metric fluctuation 
dgi' = g^'^Sguit, y, x, r)Y{xy /^i "^i ^i 0); which in the Einstein frame obeys the equation of 
motion 



D 



Sgi 



0. 



wd,Eogi = u. (E.6) 

The other Einstein equations impose additional constraints, namely Sgi must be inde- 
pendent of Xi,X2 coordinates. Substituting the spherical harmonics and the black hole 
metric into flE.6p leads to 



□, 



5d,KK 



/t(r)-3 /(/ + 4) + (A;(r)-l)n: 



Sgi 



0, 



(E.7) 



where the D'Alembert operator n^d.KK is written with respect to the appropriately 
rescaled non-compact space metric, such that the 5d Einstein-Hilbert action is prop- 
erly normalized upon performing a standard Kaluza-Klein reduction on the black hole 
background: 



d^^ 



k{r) 



k{r) 



{dt + dyf-l^dt^+ "^^ ' -^-^2 



k{r) 



dl_ 
k{r] 



+ dx' 



dr 



2 1 



r^f{r) 



. (E.8) 



(E.4) 
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In momentum space, and in Minkwoski signature the equation of motion for the metric 
fluctuation reads 

5gi'^ = F{uj,py,p = 0;u), 

F = u^{l-u)^{l + u)^H{u), 

H" +{- + + —^ ) H' + — — '\,, ^, ^, H = 0, 



u u—l u + lj u{u — l){u + l) 

1(1 + 4:) 

4 '■ 



X = l-\ l + A^iPy + Ljy 



A-/3r+, ^^^, Py^^, (E.9) 



where the i^-equation is of Heun type, with parameters 



7 = 1 - 2^/ 1 + A2(p^ + a;)2 + ^^^±^1, 6 = 1- iuje. £ = 1 + u;. 



a = P = ^— ^ + 1-^1 + K^{py + oj 



l-^. , . .A , .,.„ , ,.^2^i([±i) 



l + i 



2Wl + A2(p, + u.)2 + i(i±il-l 



^'-P«-^^- (E.IO) 



(We hope that the parameter of the Heun function, /3, is not confused with the parameter 
/3 of the black hole metric, which is now appearing through A.) 

In (lE.gp we have retained only the solution which behaves like an incoming wave at the 
horizon. The case of the trivial spherical harmonic (/ = 0) was discussed in [TT| WI\ . and 
was the template for the minimally coupled massless scalar fluctuation equation of motion 
(of course, the metric fluctuations transform as a rank 2 symmetric tensor, and this is the 
reason for the additional constraints which were placed on 5gi ). Here we are interested 
in the equation of motion obeyed by massive fluctuations in the black hole background, 
which is why we considered the most general type of scalar spherical harmonic, leading 
to the equation of motion flE.7p . Since the S^ factor is squashed in the lOd black hole 
geometry, both / and n^ quantum numbers appear in ( IE.7p . However, it is reassuring that 
despite the complexity of the metric and of the peculiarities of the equation of motion for 
the massive fluctuations, the solution can still be obtained through Heun's functions. 
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